Mon. Not. R. Astron. Soc. 000, [l]{T9] (2011) Printed 15 September 2011 (MN style file v2.2) 



HR5907: Discovery of the most rapidly rotating magnetic 
B-type star by the MiMeS Collaboration^f 

J.H. Gmnhuti'4, Th. Rivinius^, G.A. Wade^, R.H.D. Townsend^ W.L.F. Marcolino^ 
^ D.A. Bohlender^ Th. Szeifert^ V. Petit^ J.M. Matthews^ J.F. Rowe^ 
O A.F.J. Moffat^o, T. Kallinger8'i\ R. Kuschnig^'^i, D.B. Guenther^^ 

S.M. Rucinskii3, D. Sasselovl^ W.W. Weiss^^, and the MiMeS Collaboration 

iT^ ^ Dept. of Physics, Engineering Physics & Astronomy, Queens University, Kingston, Ontario, K7L 3N6, Canada 
_ ^ Dept. of Physics, Royal Military College of Canada, P.O. Box 17000, Station Forces, Kingston, Ontario, K7K 7B4, Canada 
^ESO - European Organisation for Astronomical Research in the Southern Hemisphere, Casilla 19001, Santiago 19, Chile 
"^Dept. of Astronomy, University of Wisconsin-Madison, 2535 Sterling Hall, 475 N Charter Street, Madison, WI, 53706, USA 
^— H ^ Universidade Federal do Rio de Janeiro, Observatorio do Valongo Ladeira Pedro Antonio, 43, CEP 20080-090, Rio de Janeiro, Brasil 
^National Research Council of Canada, Herzberg Institue of Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada 

Dept. of Geology & Astronomy, West Chester University, West Chester, PA, 19383, USA 
^ Dept. of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC, V6T IZl, Canada 
J/^ ^NASA Ames Research Center, Moffett Field, CA 94035, USA 

• ^^Dept. de Physique, Universite de Montreal, CP. 6128, Succursale: Centre-Ville, Montreal, QC H3C 3J7, Canada 
Institut fiir Astronomic, Universitdt Wien, Tiirkenschanzstrasse 17, A-1180 Wien, Austria 



^'^ Institute for Computational Astrophysics, Dept. of Astronomy and Physics, Saint Marys University, Halifax, NS, B3H 3C3, 



^ ^'^Dept. of Astronomy and Astrophysics, University of Toronto, 50 St Ceorge Street, Toronto, ON M5S 3H4, Canada 
» I Harvard— Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 

C3 



Canada 



15 September 2011 



> 

in 



ABSTRACT 

We report the discovery and analysis of a very strong magnetic field in the rapidly 
rotating early B-type star HR5907, based on observations obtained as part of 
the Magnetism in Massive Stars (MiMeS) project. We infer a rotation period of 
0.508276l;Q QgQg}2 d from photometric and Ha EW measurements, making this the 
shortest period, non-degenerate, magnetic massive star known to date. From the com- 
I parison of lUE UV and optical spectroscopy with LTE bruce/kylie models we find 

a solid-angle integrated, uniform black-body temperature of 17000 ± 1000 K, a pro- 
.5^ jected rotational velocity of 290 ± lOkms"^, an equatorial radius of 3.1 ± O.2i?0, 

a stellar mass of 5.5 ± 0.5 M0, and an inclination angle of the rotation axis to our 
^ line-of-sight of 70 ± 10° . Our measurements of the longitudinal magnetic field, which 

vary between -500 and -2000 G, phase coherently with the rotation period and imply 
a surface dipole field strength of '--^15.7kG. On the other hand, from fits to mean 
Least-Squares Deconvolved Stokes V line profiles we infer a dipole field strength of 
'^10.4kG. This disagreement may result from a magnetic configuration more complex 
than our model, and/or from the non-uniform helium surface abundance distribution. 
In either case we obtain a magnetic obliquity nearly aligned with the rotation axis 
{j3 = 7^^°). Our optical spectroscopy also shows weak variability in carbon, silicon 
and nitrogen lines. The emission variability in hydrogen Balmer and Paschen lines 
indicates the presence of a dense, highly structured magnetosphere, interpreted as a 
centrifugally supported, magnetically confined circumstellar disk. 

Key words: stars: individual HR5907, stars: magnetic fields, stars: rotation, stars: 
circumstellar matter, techniques: photometric, techniques: polarimetric 



* Based on observations obtained at the Canada-France-Hawaii 



Telescope (CFHT) which is operated by the National Research 
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1 INTRODUCTION 



2 OBSERVATIONS 



Since stars more massive than ~2 AIq lack a significant con- 
vective envelope, which is a necessary ingredient to drive a 
solar-type dynamo, magnetic fields are not expected to be 
observed in these stars. However, the chemically peculiar, 
intermediate mass Ap/Bp stars have been known for over 
half a century to host strong, globally organised magnetic 
fields with surface strengths of up to tens of kG (e.g. |Borra 



fc Landstreet|[l973 K At the hotter end of the Bp star class 



are the main sequence. He-strong stars (e.g. Bohlender et al 



1987 1 that show significant enhancement and often variabil- 



ity in their helium lines, as is found in the archetypical star 
a Ori E (e.g. [Landstreet fc Borra||1978[) an d the recently 
discovered star HR7355 ( Oksala et al.|2010 Rivinius et al 
2010). In addition to the helium variability, some He-strong 



stars show emission variability in Balmer lines, photomet- 
ric brightness variations, variable UV resonance lines, and 
non-thermal radio emission, most of which vary with a sin- 
gle period, interpreted to be the rotational period of the 
star (e.g. jPedersen fc Thomsen||1977[ |Walborn|1982[ [Shore 



fc Brown|fl990| [Leone fc Umana||1993i . Many of these phe 
nomena are thought to be due to the presence of a rigidly 
rotating, centrifugally supported magnetosphere - a region 
in the circumstellar environment where the stellar wind cou- 
ples to the magnetic field and is forced to co-rotate with the 



star ( [Shore fc Brown|1990[ [Townsend fc Qwocki[[2005 1 

The subject of this paper, HR5907 (H P 142184, V1040 
S co), is a bright ( V=5.4), early type B2.5V ([Hoffleit fc WaF 



ren|l991 1, emission line star with a high v sin i (340 km , 
Fremat et al. 20051. HR5907 is located in the nearby 
Upper Scorpius OB association at a distance of ~145pc 
( [Hernandez et al.|2005t . 

This star was originally brought to the attention of the 



Magnetism in Massive Stars (MiMeS) collaboration (Wade 
et al. 20111 because of its P ~ 0.5 d Hipparcos photomet- 
ric variability ( [Hubert fc Floquet[|1998| ), which would be 
in good agreement with the expected rotation period for a 
star of this spectral type and «e<j sini. Additionally, archival 
FEROS spectra suggested that the observed Ha profile is 
morphologically more similar to the profiles of other mag- 
netic He-strong stars (such as a Ori E or HR 7355) than to 
classical Be stars. 

In this paper we report on the discovery and first de- 
tailed investigation of the magnetic field and spectral vari- 
ability properties of HR5907. In Sect. [2| we discuss the de- 
tails of the observed spectroscopic, spectropolarimetric and 
photometric data, while in Sect.|3[we re-examine the photo- 
metric and Ha variability to refine the period. The funda- 
mental parameters are evaluated in Sect. [4] and the spectral 
variability is analysed in Sect.[5| In Sect. [6[ we constrain the 
magnetic field strength and geometry using the polarimetric 
data. In Sect. [7| we examine the properties of the magneto- 
sphere and present a discussion of this study in Sect. [8| 



Council of Canada, the Institut National des Sciences de I'Univers 
of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. 

t Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, 
under Prog-ID 284.D-5058. 
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2.1 Polarimetry 

Between 2010 February and 2011 March, 27 high-resolution 
{R ~ 68 000) spectropolarimetric (Stokes / and V) ob- 
servations of HR5907 were collected with the ESPaDOnS 
spectropolarimeter at the Canada-France-Hawaii Telescope 
(CFHT) as part of the Survey Component of the MiMeS 
Large Program. Each spectropolarimetric observation con- 
sists of four individual sub-exposures that were processed 
using the upena pipeline running libre-esprit, following 
the double-ratio procedure as described by |Donati et al.| 
( 1997 1 to produce the Stokes I and V spectra. Null polari- 



sation spectra were also produced by combining the individ- 
ual four sub-exposures in such a way that the polarisation 
should cancel out ( [Donati et al. 19971. This allows us to 
verify whether spurious signals are present in our reduced 
data. 

Low-resolution {R ~ 5 000) spectropolarimetric ob- 
servations were also collected with the FORS2 spec- 
tropolarimeter at the European Southern Observatories' 
(ESO) Very Large Telescope (VLT). Over six nights 
between February and May 2010, forty-eight exposures 
were collected, eight exposures per night, yielding six in- 
dividual observations. The sequence of eight exposures 
were taken with the retarder plate positioned at angles 
of -45°, -45°, -H45°, -h45°, -45°, -45°, -^45°, +45° with re- 
spect to the axis of the WoUaston prism. The Stokes V 
parameter was extracted for each observation following the 



method of [Bagnulo et al. ( 2002 1 . A log of all our polarimetric 
observations is given in Table [jj 

To increase the signal-to-noise ratio (SNR) of our ES- 
PaDOnS spectra, we applied the Least-Squares Deconvo- 



lution (LSD) procedure of [Donati et aL\ ( [1997[ ) to all ES- 
PaDOnS polarimetric data. As a starting point, our mask 
was based on a simple B2.5 star template containing only 
lines with intrinsic line depths greater than 10 percent of the 
continuum. We then proceeded to remove all lines that were 
blended with hydrogen lines and any that were too weak to 
be visible in the spectrum of HR 5907 due to the high rota- 
tional broadening of the lines. What remained is a mask of 
88 lines, most of which are helium lines. The resulting LSD 
profiles were computed on a spectral grid with a velocity 
bin of 9.0kms~^, an example of which is shown in Fig. [ll 
The LSD method improved our sensitivity to weak Zeeman 
signatures by increasing the SNR by a factor of about 11 
times compared to the peak SNR. We also list the improved 
LSD SNR measurements for each observation in Table[l] We 
find that each observation shows a clear Zeeman signature 
in the Stokes V profile and all but the profiles obtained on 
2011-03-16 have a definite detection (false al arm probabil- 
ity (FAP) < 10~^) according to the criteria of Donati et al 



( 2002[[2006[ ). The two profiles from this night have FAP of 
1.677 X 10"^ and 7.744 x 10"^ In no case do we detect a 
significant signal in the null profile for any of our extracted 
mean profiles (FAP > 10"^). 

The longitudinal magnetic field (Be) and null measure- 
ments from the ESPaDOnS spectra were computed from 
each LSD Stokes V and diagnostic null profile in a manner 



similar to that described by Silvester et al. (20091, using 



an integration range from -400 to 400 km s . We find that 
the Be varies between -740 and -1845 G with a typical un- 
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Table 1. Journal of polarimetric observations listing the date, the heliocentric Julian date (2,455,000+), the number of sub-exposures and 
the exposure time per individual sub-exposure, the phase according to Eq.[l] the peak signal-to-noise ratio (SNR) (per 1.8 and 3.4kms~^ 
velocity bin in the ESPaDOnS spectra and FORS spectra, respectively) in the unpolarised spectra, the mean SNR per 9.0kms~^ velocity 
bin in the LSD Stokes V profile, and the derived longitudinal field and longitudinal field detection significance z from both V and A'^. 
All ESPaDOnS observations, except the profiles obtained on 2011-03-16, have a definite detection of a signal in Stokes V (FAP < 10~^) 
while no detection of a signal in the diagnostic null profile is found (FAP > 10"'^). 



HJD tcxp Peak LSD V N 

Date (2455000-I-) (s) Phase SNR SNR ± crs z Bt±aB z 
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certainty of ~120G. The longitudinal magnetic field was 
measured from the helium lines in the FORS data and were 



obtained in a manne r similar to Rivinius et al. (20101 and 
Bagnulo et al. (2009 Appendix 2) using a minimisation 



technique. The uncertainties were computed following the 



method outlined by Rivinius et al. ( 2010 1 using a bootstrap 



Monte Carlo approach. We find a very similar range of values 
from -750 to -2048 G, but with typically larger uncertainties 
of ~250G. The longitudinal field measurements from each 
polarimetric observation are also listed in Table [l] 



2.2 Spectroscopy 

A number of high-resolution spectra were acquired with the 
Ultraviolet and Visual Echelle Spectrograph (UVES) at the 
VLT. The instrument was used in its DIC2 437/760 setting 
giving us blue spectra from 375 to 498-nm and nearly con- 
tinuous red spectra from 570 to 950-nm. A slit width of 0.8" 



was used resulting in a resolving power of = 50 000 over 
the entire spectrum. The details of these observations are 
listed in Tabled 

In addition to the observations already described, 
archival lUE UV data obtained via the NASA-MAST 
archive were also used to constrain the fundamental pa- 
rameters of HR5907 (see Sect.|4|. High-resolution and low- 
resolution observations exist in both the SWP and LWP 
spectral regions, both taken in large aperture mode. The 
high-resolution and low-resolution spectra are in good agree- 
ment, but we have only used the low-resolution spectra (both 
obtained on the night of 1989-07-26) here due to the higher 
SNR compared to the high resolution spectra. 



2.3 Photometry 

Between 2011 April 15 and 2011 May 03 we obtained ap- 
proximately uniform sampling of the photometric brightness 
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Table 2. Journal of UVES observations listing the date, the heliocentric Julian date (2,455,000+), the exposure time, the phase according 
to Eq. rn and the peak signal-to-noise ratio (SNR) per 2.4kms~^ velocity bin. 



HJD texp Peak HJD toxp Peak 

Date (2455000+) (s) Phase SNR Date (2455000+) (s) Phase SNR 
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Figure 1. Mean LSD Stokes V (top), diagnostic null (middle) 
and Stokes / profiles (bottom) of HR 5907 from 2010 February 24. 
The V and N profiles are expanded by the indicated factor and 
shifted upwards for display purposes. A clear Zeeman signature is 
detected in the Stokes V profiles, while the null profile shows no 
signal. The error bars shown in the Stokes V and diagnostic null 
profiles represent the la uncertainties. The integration limits used 
to measure the longitudinal field in the ESPaDOnS LSD profiles 
are indicated by the dotted lines. 



variations of HR 5907 using the Microvariability and Oscil- 
lations of STars (MOST) satellite. The observations were 
obtained with the MOST satellite in switched-target mode, 



meaning that observations were scheduled such that mul- 
tiple targets were observed during each MOST orbital pe- 
riod. During our observing run, half of the orbital phase was 
shared with observations of Arcturus. Each observation of 
HR5907 consisted of 157 stacked 0.18 s exposures, resulting 
in a total of 9846 observations. 

Photometric measurements of HR5907 were obtained 
by performing aperture photometry on 20 x 20 pixel sub- 
rasters obtained from the MOST CCD photometer. We used 
a radius of 3 pixels for the photometric radius and a sky 
annulus of 8 pixels. Centroids were measured by fitting a 
Gaussian Point-Spread-Function (PSF). As with aU MOST 
photometric studies, the raw photometry shows strong pho- 
tometric variations correlated to the 101 minute MOST or- 
bital period. The variations are attributed to an increase 
in the measured background level of ~200 ADU and to 
centroid shifts that map out intrapixel sensitivity changes 
across a CCD pixel. The photometric variations attributed 
to changes in the background level were measured to be 
~2 percent and changes in centroid position to be 0.3 per- 
cent. HR5907 also shows intrinsic periodic variations of ~4 
percent. In order to extract photometry free of instrumental 
effects a simultaneous fit was performed. The stray fight was 
modelled with a third order polynomial, the intrapixel varia- 
tions were modelled as a linear trend and intrinsic variations 
were modelled as a sinusoidal function with five additional 
harmonics. The fit was iterated after removing values of the 
centroid position and the FWHM of the PSF beyond the 90 
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Figure 2. Brightness variations of HR5907 as measured by tlie 
MOST satellite (black circles). The symbol size is chosen to rep- 
resent the mean uncertainty of the data. The inset shows an ex- 
pansion of the MOST data over a few nights to highlight the cycle 
to cycle variability and the clear periodicity. Also included is a 
red curve that corresponds to a best fit to the phased data (see 
Sect.|3]for further details). 



percentile. All data was excluded when the satellite passed 
through the South- Atlantic- Anomaly due to the large num- 
ber of cosmic ray hits on the detector. The final reduced 
data set had an effective duty cycle of 26 percent with 6062 
photometric points. 

The 18 days of observations cover roughly 36 cycles of 
the star, as illustrated in Fig. [2] These observations clearly 
demonstrate the approximately 0.5 day rotational cycle (see 
Sect. |3] and Fig. [s] for further details). Also shown in Fig. [5] 
is an expansion of a few cycles that emphasises the non- 
sinusoidal, but periodic brightness variations. The sharpness 
of the photometric minimum is not consistent with bright- 
ness patches or other surface features and requires a struc- 
ture that is geometrically thin, likely of circumstellar origin. 
We note that the observed photometric variations are consis- 
tent from cycle to cycle, but that there are additional small 
intrinsic variations beyond the expected scatter about the 
uncertainty, which cannot be attributed to any instrumen- 
tal or orbital effects. 



3 EPHEMERIS 

HR5907 was first identified as a photometric variable star 
with a period of 0.508 d based on Hipparcos data ( |Ifubert"fc| 
Floquet 1998 1. Lefevre et al. ( 2009 1 later confirmed this vari- 



ability but adopted a period that was twice as long (1.0 17 d), 
since they assumed it was due to binary effects. We re- 
analysed the Hipparcos dataset using a Lomb-Scargle like 
technique ( Press et al.|1992 I and confirm a best-fit period of 
0.50831 ± 0.00003 d. We also confirm that the current epoch 
of longitudinal magnetic field measurements. Ho equivalent 
width (EW) variations (see Sect. [5| and MOST photom- 
etry phase well with this period. However, because of the 
long baseline between the Hipparcos era photometry and 
the current MOST photometry, choosing any period within 
the uncertainty range can lead to a phase shift between the 
phased Hipparcos and MOST photometry of upwards of 0.5 



cycles. Unfortunately, as the MOST data only covers about 
36 rotation cycles, we cannot obtain a more precise period 
from this dataset. Applying the Lomb-Scargle like technique 
to our extensive Ha EW dataset does provide a more precise 
period of 0.50825 ±0.00001 d, but this period is inconsistent 
with the Hipparcos period. In either case, we believe that 
the periodic variations are sufficiently non-sinusoidal that 
the Lomb-Scargle technique is inadequate. 

We proceeded by adopting a multi-harmonic fitting 



technique similar to Schwarzenberg-Czerny ( 19961, which is 
ideal for non-sinusoidal periodic variations. After including 
contributions from the first three harmonics to both the Ha 
EW variations and the Hipparcos photometry we find a best- 
fit period of 0.508274l|5 E;j5°oi2 d fro™ the EW varations and 
a best-fit period of 0.508276lE;:Si2 d from the Hipparcos 
photometry. A similar analysis of the MOST data results 
in a best-fit period of 0.508269l^ J5E!J!!54j d period. The uncer- 
tainties on these newly derived periods are still large enough 
that the phasing uncertainty between the Hipparcos-era and 
current epoch of data is not resolved. In any event, we adopt 
the Hipparcos period as it provides good agreement between 
the Hipparcos and MOST photometry, as is expected unless 
significant rotational braking has occurred (see Sect, [s] for 
further details). Therefore, within the context of an oblique 
rotator model, we adopt this period as the rotational period 
of HR5907, which results in this star having the shortest 
known rotational period of any non-degenerate, magnetic 
massive star. Taking the Ha EW maximum as To we derive 
the following ephemeris: 

HJDS'^" = 2447913.694(1) + 0.508276(tl^) ■ E, (1) 

where the uncertainties in the last digits are indicated in 
brackets. Unless otherwise stated, all further data are phased 
according to this ephemeris. 



4 FUNDAMENTAL PARAMETERS 

In order to investigate the fundamental parameters of 
HR5907 we utilised the third revision of the BRUCE 
and KYLIE software suite (hereafter referred to as BK3; 



Townsend 19971. While BK3 is capable of producing the- 



oretical spectra in both absolute flux and continuum nor- 
malised flux from either ATLAS9 LTE atmo spheres (Kurucz 
[T992) or from NLTE tlusty atmospheres of Lanz & Hubeny 
^2007j) , only solar abundance LTE atmospheres were used in 
this study. This was decided due to the fact that we found 
a much better overall agreement between the line strengths 
of C 11 and Si ll lines in the blue region with solar abun- 
dance LTE atmospheres in comparison to solar abundance 
NLTE atmospheres. In our analysis, both the observed spec- 
tra and synthetic spectra were renormalised to the contin- 
uum regions surrounding the individual line profiles before 
the comparisons. 

In the following discussion, we assume that modelling a 
rapidly rotating star depends on five independent parame- 
ters, which we choose to be the equatorial velocity Veq sini, 
the inclination angle of the rotation axis relative to the line- 
of-sight i, the effective temperature Tofi, the stellar mass 
Mi,, and the stellar equatorial radius -Reg. As input, BK3 
uses our estimated period P = 0.508276 d to constrain the 
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Figure 3. Phased observational data according to ephemeris given in Eq. [T] Upper panel: MOST (black triangles) and Hipparcos 
(red squares) photometry. The MOST data have been binned in 0.01 phase bins and shifted to the Hipparcos magnitude for display 
purposes. Middle panel: longitudinal magnetic field variations measured from ESPaDOnS (black circles) and FORS (red diamonds) 
data. Lower panel: Ha equivalent width variations measured from ESPaDOnS (black circles) and UVES (green triangles). We have 
also drawn dotted lines indicating phases 0.0 and 0.5, and dashed lines indicating phases 0.25 and 0.75. 



equatorial values for relevant parameters. In order to de- 
termine the properties of HR5907, we searched through a 
grid of models that varied the previously listed fundamental 
parameters for a fixed P. 



4.1 Effective Temperature 

Our primary method to constrain the Tcfi of HR5907 was 
from a de-reddened Spectral energy distribution (SED). We 
note that BK3 uses T^a in the sense that it is the solid angle 
integrated, uniform black-body temperature that a star of 
the same surface area would need to have, such that the total 
luminosity is the same as the actual gravity darkened star. 
In this context, BK3 allows for different local temperatures 
over the surface of the star. 

We constructed an SED for HR5907 by combining 
the lUE UV data with flux calibrated UBV photometry 



The SED was then extinction corrected using the para- 
metric law of Cardelli, Clayton & Mathis (19891, with 



( Jaschek fc Egret||1982| [Kharchenko fc Roeser|[2009t , using 
calibrations of |Hynes| ( |2011[ ). To de-redden our spectra we 
adopted a {B — V)o = —0.18 corresponding to a temperature 
between 17 to 18 kK. This temperature, as suggested by our 
initial fits to the Cll 4267 A line, is cooler than otherwise 
suggested by IIR5907's B2.5 spectral classification. How- 
ever, based on our initial temperature fits and the stronger 
than predicted helium absorption lines, HR 5907 may in fact 
be a He-strong star and therefore the cooler temperature is 
justified. With this [B - V)o we find that E{B -V) = 0.14, 
which is slightly lower than the value determined by |Pa-| 
paj, Wegner fc KrelowskI] (p9T|) {E{B - V) = 0.155) who 



adopted a standard B2 star template for their analysis. We 
do not find any significant differences in our results if we 
adopt this slightly higher value. 



Rv ~ 3.1, and we then normalised the SED by the fiux 
at 5500 A so that we could constrain the temperature by fit- 
ting the slope of the UV spectrum, removing the need for any 
distance or luminosity corrections. As illustrated in the top 
panel of Fig.|4] we find a best-fit T^s = 17000 ± 1000 K when 
fitting the 1250-2000 A range of the SED. We note that if we 



adopt an Rv = 3.9 as found by Lewis et al. (20091, we find 



a slightly cooler temperature of 16500 K, but still within our 
uncertainty. We also confirm that adopting a NLTE model 
atmosphere results in a best-fit T^s within our uncertainty. 



4.2 Projected Rotational Velocity and Radius 

As previously mentioned, we could also independently con- 
strain the Tea by fitting the hue depth of the Cll 4267 A 
line profile, as i or has very little affect on the shape 
of this line and it only varies slightly with R^q. However, 
there is still a strong dependence of the line depth from the 
rotational broadening. 

We simultaneously fit the T^s and the projected ro- 
tational velocity Veq sin i of the star by searching for the 
model that provided the best overall fit to the C 11 profile. We 
first constructed a mean profile from all of our ESPaDOnS 
spectra before comparing this profile to our grid of mod- 
els. The ESPaDOnS profile was best fit by a model with 
Tefi = 17000 ± 1000 K and v^q sini = 290 ± 10 km s~^ 

As a check, we also compared our models with a mean 
profile constructed from the UVES data that resulted in 
a best-fit Tcfi = 17000 K, and a slightly lower Ue^sini = 
285kms~^, which is still consistent with our findings from 
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Figure 4. Comparison between BK3 model and observed data: Upper panel: lUE UV spectrum normalised to flux at 5500 A and 
best-fit model with T^g of 17000 K (red) along with models corresponding to 16000/18000 K (blue/green). Middle panel: lUE UV 



an inclination of i = 70°. Bottom panel: ESPaDOnS normalised spectrum (corresponding to observation with the least emission (24 
February 2010; black)) with best-fit overall model (red). Note that the poor fit to the Balmer wings (bottom panel) is likely due to the 
over-normalisation of the ESPaDOnS spectrum. 
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Figure 5. Comparison between the mean profile from ES- 
PaDOnS data (black) and best-fit model with Teff = 17000 K, 
vs'mi = 290kms~^ (red) and vsini = 280 and 300kms~^ 
(dashed blue and green, respectively). 



the ESPaDOnS spectra. We adopt Veq sini = 290±10kms~^ 
since our mean ESPaDOnS profile is constructed from a 
set of observations that better span the entire rotational 
cycle. A comparison between our best-fit model and our 
mean ESPaDOnS profile is shown in Fig. [S] By combining 
our best-fit Veq sin i and rotational period P, we find that 



,sini ~ 2.91 i?, 



,0. 



4.3 Stellar Mass and Inclination 

With Veq sin i and TcH determined, the remaining unknown 
model parameters are M*, and i, since Rsq is constrained 
by P and VeqSini. To infer the mass of HR5907 we made 



(tt — 7.64 ± 0.37; van Leeuwen 20071 to obtain the abso- 
lute magnitude My and therefore obtain HR5907's place- 
ment on a Hertzsprung-Russell (EtR) diagram. The visual 
extinction Ay and magnitude V were taken to be 0.43 (this 



work) and 5.4 (Kharchenko & Roeser 2009 1, respectively. Us- 
ing a bolometric correction corresponding to Teg from our 



model {BC = -1.6 ± 0.2; |Lanz fc Hu"beiyl[2007l ) we com- 
puted the bolometric magnitude (Mbol ~ My + BC) and 
finally the luminosity L*. Our final luminosity is found to 
be log(Z/*/L0) = 2.78 ±0.14 (using the absolute bolometric 
corrected solar magnitude Mq — 4.74), with the uncertainty 
found through propagation of the uncertainty in the distance 
and the bolometric correction. 

As shown in Fig. |6] we compared HR 5907's placement 
with available Padova evolutionary tracks ( y=0.26; Bertelli 
et aH2009) and find that HR 5907 has a mass of ^ 5.2 Mq. 
Also included in Fig. |6] are the results of ( Hernandez et al. 



2005 1 who find HR 5907 slightly hotter and more luminous, 



but still consistent with our findings. The Padova evolution- 
ary tracks may lack any effects due to rotation, but accord- 
ing to the work of^Maeder & Meynet ( 2010 1, the inclusion of 



rotation would shift the evolutionary track such that a star 
of a particular mass and age would be slightly cooler and 
slightly brighter when rotation is included. Therefore, if we 
used tracks that included rotation, then based on HR 5907's 
position we would only expect to find small changes (within 
our uncertainty) that would result in this star being slightly 
closer to the main sequence and having a slightly higher 
mass. 

We can also constrain the mass and inclination from fits 
to the Balmer wings, which are sensitive to surface gravity 
{g = G Mi, / R^). Unfortunately, this can be quite problem- 
atic with HR5907 as there is considerable emission in the 
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Figure 6. Left panel: Hertzsprung-RusscU diagram showing tlie 
position of HR 5907 as found in this work (black square) and as 
found by |Hernandez et al.| | |2005l l (red circle). Also shown are stel- 
lar evolutionary tracks (dashed) and the zero-age main sequence 
(solid) from [Bertelli et al.| l |2009l l. Right panel: Comparison of 
our best-fit model (T^s = 17kK, M* = 5.5Mq, Req = 3.1i?Q) 
with the FORS observation from 25 April 2011 showing the least 
emission, for the H(5 region. 



wings of the Balmer lines due to circumstellar plasma (see 
Sect. [tJi. However, we still proceeded by using the FORS 
spectrum from 25 April 2011 that showed the least emission 
in the Balmer wings. The long-slit FORS spectra are ideal 
since they are less susceptible to the over-normalisation of 
the Balmer lines that can occur with the multi-order ES- 
PaDOnS and UVES spectra. We proceeded by searching for 
the model (at a fixed Toh , but varying mass and radius) that 
provided the best fit to the wings of the H(5 line profile, since 
lower Balmer series lines are more affected by the circumstel- 
lar plasma, and higher Balmer series lines were more difficult 
to ensure consistent normalisation between the model and 
observed spectrum. An overall best-fit model was found with 
Mi, = 5.5 ± 0.5 Mq and R^q = 3.1 ± 0.1 Rq (a comparison is 
shown in the right panel of Fig.|6|, consistent with the mass 
estimates derived from HR 5907's placement on the HR dia- 
gram. If we adopt = 5.5±0.5 and Req = S.liO.li?© 
this implies that values of i = 70tio° allowed. 

Independently, we can also use the BK3 models to fit the 
absolute, extinction-corrected SED of HR5907. Using the 
previously stated parallax value, we first corrected the SED 
to an absolute distance consistent with the BK3 models. 
We then proceeded by fitting the lUE UV spectra (between 
1250-2000 A) with the BK3 models. From this method we 



also find a best-fit R^q = 3.1 ± 0.1 Rq and i = 70 ± 10° for 
masses around 5 Mq , as shown in the middle panel of Fig. [1] 
However, if we use the Rv value of Lewis et al. ( 2009[ ) and 
the slightly cooler temperature that results from those fits, 
we would find a slightly smaller equatorial radius of 2.9 i?0. 
We therefore adopt a slightly higher uncertainty of Q.2Rq. 
A summary of the fundamental stellar parameters is listed 
m Sect. [8] A comparison between our best-fit model and the 
ESPaDOnS spectrum with the lowest emission is included 
in the bottom panel of Fig. [4] 



5 LINE VARIABILITY 

Because of HR 5907's rapid rotation and the need to sample 
fast photospheric variations, we opted to use the individual 
sub-exposures that make up the polarimetric ESPaDOnS 
observations for all our line profile variability analysis. The 



individual sub-exposures are still of sufficient quality, with 
SNRs comparable to the UVES data. 

To begin, we can characterise the variability from 
variations in the equivalent widths (EWs) of the spectral 
lines. Before measuring the EW, each spectral line was re- 
normalised to the surrounding continuum, and a telluric cor- 
rection algorithm was applied to regions redward of 5790 A. 
The EWs were then obtained by numerically integrating 
over the line profile. The la uncertainties were computed 
by adding the individual pixel uncertainties in quadrature. 
For the UVES data, a single uncertainty value was assigned 
to each pixel, determined from the RMS scatter in the con- 
tinuum regions surrounding the line profile. The EW varia- 
tions for a number of spectral lines were measured and then 
phased with the ephemeris given in Eq. [l] with the results 
presented in Figs. |3] [T] and |8] We find small systematic off- 
sets between the UVES data and the ESPaDOnS data on the 
order 0.01 A for most helium and metallic lines. However, a 
larger offset is found in the broader Balmer and Pashen lines. 
We do not find any significant differences in the morphology 
of the line profiles between the two datasets at the same 
phase, but we do find weak but wide-scale systematic differ- 
ences. We attribute these systematic offsets to inconsisten- 
cies in the normalisation between the two datasets, which is 
common due to complexity of extraction and normalisation 
of echelle spectra. Therefore, the UVES EW values have all 
been corrected to best match the ESPaDOnS measurements 
by subtracting fixed offsets between the two datasets. 

As depicted in Fig. [7j helium lines show obvious EW 
variations, consistent with the rotational period, as mea- 
sured from our high-resolution spectra. The observed EW 
variations do not vary in a similar manner to the circumstel- 
lar Hq variations depicted in Fig. [S] and we speculate that 
they are due to photospheric abundance patches, consistent 
with findings from investigations of other He-strong/variable 
stars (e.g. Bohlender fc Landstreet|1988 1. 

The clearest variability is seen in the He i 4920 A hue, 
which reaches maximum absorption (corresponding to a 
minimum EW in Fig. [7| around phase 0.25 and minimum 
absorption around phase 0.55. Another local minimum is 
found around phase 0.85, indicating a relatively complex he- 
lium abundance pattern in the photosphere of HR 5907. In 
contrast to the large EW variations seen in He I 4920 A, very 
little variability is observed in the forbidden [He l] 4045 A 
line. In the other helium lines shown in Fig. [7] we see a vari- 
ability pattern similar to the 4920 A line, but with a smaller 
amplitude. In comparison to other He-strong stars such as 
0- Ori E (e.g. [Landstreet fc Borra|1978|[Reiners et al. l2000| 



Smith & Bohlcndcr''2007|) or HR7355 flRivinius et al.|2010| 



we find that HR 5907 shows relatively weak variability in its 
helium lines, of ~20 percent in EW. 

We also investigated the EW variability of a number 
of other photospheric absorption lines, the results of which 
are displayed in Fig. |8] In comparison to the helium lines, 
very little variability is seen in these lines. Only the lines of 
carbon and nitrogen, as evidenced in Fig. [s] (top left and top 
right panels) show any obvious signs of periodic variations. 
The EW curve for carbon is approximately sinusoidal, with 
maximum absorption around phase 0.2 and minimum ab- 
sorption at about phase 0.7. The nitrogen EW curve is also 
simple, but appears to be a non-sinusoidal variation. 

We can further investigate the photospheric variations 
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Figure 7. Phased equivalent width variations of selected helium lines in the spectrum of HR 5907, measured from ESPaDOnS (black 
circles) and UVES (red triangles) data. 
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Figure 8. Phased equivalent width variations of selected metal lines in the spectrum of HR 5907, measured from ESPaDOnS (black) 
and UVES (red) data. 



by examining the phased intensity variations of these line 
profiles, commonly referred to as "dynamic spectra" , as 
shown in Figs. [9]and |10[ In these figures, individual spectra 
are represented as horizontal bands with intensity encoded 
by different colours and are stacked vertically according to 
their phase to illustrate time variability. To further highlight 
these variations, we have also subtracted from each profile 
the mean spectrum formed by averaging all the profiles. 

The dynamic spectrum of the [He l] 4045 A line (Fig. [9] 
upper left panel) shows a very weak but simple pattern 
with one clear enhanced absorption feature travelling from 
negative to positive velocities, crossing Vsys = Okms~^ at 
about phase 0.75 (where Vsys represents the systemic ve- 



locity of HR5907). The absorption feature is likely a re- 
gion of the stellar photosphere where helium is found to be 
over-abundant compared to the mean helium surface distri- 
bution. This feature is also present in all the other helium 
lines shown in Fig. |9j but appears to be the less promi- 
nent absorption feature in the other helium lines. A signifi- 
cantly larger and more absorptive feature is found to cross 
Vsys ~ 0kms~^ at about phase 0.25 in all the other helium 
lines. Additionally, these stronger helium lines also show 
a clearly more complex dynamic spectrum. In these lines 
there also appears to be the presence of a strong pseudo- 
emission feature (this feature only appears in emission rela- 
tive to the mean profile and represents a region of the stel- 
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lar photosphere where hehum is under-abundant with re- 
spect to the mean helium surface distribution), which crosses 
Vays = Okms"'^ around phase 0.6. In some of the lines (e.g. 
He I 4713, 5015, and 6678 A) additional pseudo-emission fea- 
tures can been found with one crossing Vsys = Okms~^ 
around phase 0.8 and the other around phase 0.0. These two 
features are much narrower than the feature that crosses at 
phase 0.6, and do not reach the same intensity of pseudo- 
emission. 

If we look at the dynamic spectra of the metallic lines 
(Fig. 10 1, we find a single extremely weak pseudo-emission 



feature that crosses Vsys ~ Okms~^ around phase 0.2 in Cll 
or an enhanced absorption feature crossing Vgyg — Okms~^ 
around phase 0.6 (upper left panel). We have also tentatively 
identified a weak absorption feature in N ll and Ne I crossing 
Vsys = Okms~^ around phase 0.2, which appears to be in 
anti-phase with the absorption feature found in the C ll line. 
However, these lines are very weak, making it difficult to 
distinguish any features relative to the noise level. In the 
dynamic spectrum of Sim 4553 A (bottom panel), we find 
two absorption features closely spaced in phase, with one 
crossing Vsys ~ Okrns"^ around phase 0.65 and the other 
crossing around phase 0.8. These variations appear to be in 
phase with the C ll variations. 



6 MAGNETIC FIELD GEOMETRY 

To investigate the magnetic field geometry of HR 5907 we as- 
sume that the field can be described by the dipole Oblique 
Rotator Model (ORM). This model is characterised by four 
parameters: the phase of closest approach of the magnetic 
pole to the line of sight 00, the inclination of the stellar ro- 
tation axis i, the obliquity angle between the magnetic axis 
and the rotation axis /3 and the dipole polar strength Bd- 
Our first approach was to model the longitudinal field curve, 
as shown in Fig. [s] This was carried out using a minimi- 
sation to compare the observed longitudinal field curve to 
a grid of computed longitudinal field curves to determine 
Bd and /3. For this analysis we have assumed i — 70° (as 



determined in Sect. 4.31 and a limb darkening coefficient of 
0.4. 

The resulting landscape from our fits is shown in 
Fig. [TT] Our results indicate that we are viewing the nega- 
tive (southern) hemisphere of the magnetic field and that 
HR5907 hosts a field with a polar surface strength of 
G, an obliquity angle oi — 7^2° ^nd that 



Bd = 157001^E;|] 



(f>o ~ 0.48. Despite relatively small formal uncertainties our 
best-fit model is not capable of fitting the measured longi- 
tudinal field measurements at each phase, as illustrated in 
Fig. 12 and as indicated by the best-fit reduced x^ ~ 2.0 
from the FORS and ESPaDOnS datasets combined, or 2.1 
from ESPaDOnS data alone. 

Still within the framework of the ORM, our second pro- 
cedure compares the observed mean LSD Stokes V profiles 
to a large grid of synthetic profiles that are parametrised 
by Bd, /3, the phase (ji and i. The models are computed by 
performing a disk integration of local Stokes V profiles as- 
suming the weak field approximation and a uniform surface 
abundance. The parameters of the Stokes V profiles were 
chosen to fit the mean LSD Stokes I profile created from 
the average of all our observations. We adopted the (j)o we 
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Figure 11. landscape of dipole field strength Bd versus mag- 
netic obliquity /3 permitted by the longitudinal field variation of 
HR5907, assuming j = 70°. 



obtained from our longitudinal field curve fits, but find that 
the derived value of /3 does not vary much for small changes 
in (j)o and that the derived values of Bd only vary by a few 
hundred gauss. For each Stokes V profile we found the pa- 
rameters that provided the lowest x^ for each observation. 
Using a Bayesian framework, we then combined the x^ dis- 
tributions obtained for each observation to find the maxi- 
mum likelihood model. In Fig. [13] we compare the synthetic 
profiles of our maximum likelihood model (dotted red) with 
the best-fit model for each observation (dashed blue). The 
quality of the fits are similar, showing that a single dipole 
configuration can, for most phases, roughly reproduce the 
observed mean profiles. However, as with our fits to the lon- 
gitudinal field curve, some phases are rather poorly fit by 
the ORM model. 

The maximum likelihood model was found with Bd ~ 
IO4OOI350 G and /3 = 7 ± 1°, where the uncertainties repre- 
sent the 95 percentile range. While the formal uncertain- 
ties are quite low, we note that the individual best-fits 
to each observation have best-fit parameters ranging from 
Bd ~ 4000 - 25000 G and /3 ~ - 15°. 

Therefore, modelling of the longitudinal field curve and 
Stokes V profiles indicate a strong magnetic field nearly 
aligned with the rotation axis. However, the observations 
at different phases differ quite significantly from the overall 
best-fit model. This will be further discussed in Sect.[8l 



7 THE MAGNETOSPHERE 

Within the context of the magnetosphere model, the com- 
mon picture is that the variability that is observed in hy- 
drogen lines among Bp stars results from the circumstellar 
plasma. This plasma is centrifugally supported and magnet- 
ically confined into clouds that are forced into co-rotation 



with the host star (e.g. |Shore & Brown| 1990 


Shore et 


al.||1990||Shore||1993| 


Townsend & Owocki||2005| 


ud-Doula, 


Owocki & Townsend 


2008 


1. The circumstellar gas is seen 



in emission (usually in hydrogen and a few other lines) at 
quadrature phases, and in absorption when passing in front 
of the star. 

As with the photospheric variability, we begin by first 
analysing the EW variations in the hydrogen Paschen and 
Balmer lines. As discussed in Sect.|3j HR5907 shows strong 
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Figure 9. Phased variations of photospheric helium hnes. 
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Figure 12. Longitudinal field variations of HR5907 from ES- 
PaDOnS (black circles) and FORS (red triangles) compared to 
our best-fit model curve (dashed) with = 15700 G and fj = 7° . 
The dotted lines indicated the la limits permitted by our fits. 



EW variations due to the circumstellar emission in Hq. In 
Fig. [14] the EW variations for a number of Paschen and 
Balmer lines are also depicted. The strongest variation and 



simplest pattern is that of Ha (Fig. 14 bottom left panel and 
Fig. [3] bottom panel), which shows a single emission peak at 
phase 0.5. This basic pattern is also present in the higher 
Balmer lines such as H/3 and H7 (Fig. 14 bottom panels), 



but with much more noise. However, there also appears to be 
a second, smaller absorption dip in the H/3 EW variations at 
phase 0.6 that is not present in the other Balmer EW varia- 
tions. If we now look at the Paschen lines, which probe the 
magnetosphere at a different optical depth (Fig. 



14 



top pan- 
els), we see a much more complex EW curve, composed of 
many smaller local absorption dips in addition to the single 
large emission peak. The overall structure of the EW curve 
for the Paschen lines is similar, but the amplitude of the local 
absorption variability varies significantly among the differ- 
ent lines available in the ESPaDOnS/UVES spectral range. 
We find the Pais line to show the strongest overall variabil- 
ity, but Pai4 has less scatter in the EW variations showing 
a local maximum emission at phase 0.18, followed by a local 
minimum emission at phase 0.33. Another local maximum 
in the emission is found at phase 0.69, but there also appears 
to be some additional structure to the EW curve between 
phases 0.33 and 0.69 that is poorly sampled by our data. 

We can further our understanding of the magneto- 
spheric structure by studying the dynamic spectra shown in 
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Figure 10. Phased variations of selected metallic photospheric lines. 



Figs. [15 and |16| However, unlike the dynamic spectra shown 
in Sect. 5] these dynamic spectra only include contributions 
from the emission, as a theoretical photospheric profile has 
been subtracted from each observation. In Fig. |15| we sub- 
tracted profiles corresponding to the best-fit BK3 model de- 
termined in Sect. [4] Since these models don't extend into 



the Paschen region, we found the TLUSTY model (Lanz & 
Hubeny|2007 1 that provided the best-fit to the Balmer lines 
of our BK3 model. This best-fit model (Teft = 19 kK and 
log(f7) = 4.25) was found to reproduce the wings of the BK3 
model Balmer lines well, but has a core that is less deep 
than the same BK3 model. 

From the Ha dynamic spectra we find that the main 
occultation occurs at phase 0.0, travelling redward. If we 
interpret this feature as a dense cloud in the circumstellar 
environment, then the fact that the occultation occupies a 
relatively broad velocity range during the eclipse indicates 
the cloud is azimuthally extended. Immediately following 
the broad absorption feature is a narrower feature that is 
still in emission in Ha (as highlighted in the bottom panel 
of Fig. 15 1, but is in absorption in H7 and the Paschen lines 



of Fig. |16| This feature is travelling blueward and crosses 
Vays = Okms~^ at phase 0.21. Another occultation is clearly 
visible in the dynamic spectrum of H7 moving redward and 



crossing Vsys ~ Okms~^ at phase 0.59. Two clear emission 
peaks are also visible; the first feature appears at positive 
velocities, with the centre of the feature occurring at phase 
0.36. A second blueward feature occurs at phase 0.62. We 
have included curves in Figs. [15] and [16] that indicate the 
potential orbit for a given cloud if it is in rigid rotation. 
The curves were selected such that they pass through the 
blueward-migrating occultations and reach the same veloc- 
ity extrema as the brightest emission features in the Paschen 
lines (which approximately correspond to the Kepler radius; 
see Sect, [s] for further details). From these orbital curves 
it seems likely that these two features are the same cloud 
viewed at opposite sides of the star, but the expected phase 
difference should be 0.5 and not 0.31 as measured. The mis- 
match between these curves and the migrating features seen 
in the dynamical spectra suggests that the features do not 
trace single, discrete clouds in the rigidly rotating magne- 
tosphere, but rather arise from the combined effects of mul- 
tiple, optically thick clouds distributed over a range of az- 
imuths. The observed discrepancy between the maximum 
emission of this cloud at opposite quadratures is also qual- 
itatively explained if the cloud is optically thick. Another 
redward emission feature is also visible, reaching a maxi- 
mum redward velocity at phase 0.90 in Ha and H7. This 
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Figure 13. Selected mean circularly polarised LSD Stokes V profiles (grey circles) obtained from the ESPaDOnS dataset sampling the 
rotational cycle of HR 5907. The error bars represent the la uncertainties for each pixel. Also shown are the individual best-fit model 
profiles for each phase (blue dashed) in addition to profiles corresponding to the model that provides the global maximum likelihood 
(red dotted) obtained from fits to the mean Stokes V profiles with = 10200 G and (5 = 7°. The phase and best- fit parameters are 
also indicated for each observation. 



feature appears to be related to the blueward emission fea- 
ture that is more visible in H7 at phase 0.38. If we interpret 
the emission features resulting from two distinct emission 
clouds in the circumstellar disk, it would imply that the 
clouds have a phase separation of ~0.6. 



In the dynamic spectra of the Paschen lines (Fig. 16 1 



we still see the same basic features as found in the dynamic 
spectra of the Balmer lines - two strong emission features, 
followed by a third weaker feature. However, the phases at 
which the emission features occur are not consistent with 
the phases as observed in the Balmer lines. The centre of 
the first redward emission feature is observed at phase 0.23, 
the blueward feature is found to occur at phase 0.66, and 
the last redward feature is found at phase 0.91. The differ- 
ences in the relative phasing of the features may be a result 
of subtracting the TLUSTY model and not the correspond- 
ing BK3 model, but this is unexpected as the locations of 
these emission peaks are far in the wings where the TLUSTY 
models provide the best fit to the BK3 models. There also 
appears to be a large disagreement in the phase difference 
between maximum emission, which cannot be due to using a 
TLUSTY model, as this is a constant spectrum subtracted 
from all models. The phase difference between the first red- 



ward feature and the blueward feature is 0.31 for the Balmer 
lines, but 0.43 in the Paschen lines. The phase difference be- 
tween the blueward feature and the second redward emission 
peak is only 0.25 from the Paschen lines, compared to 0.28 
from the Balmer lines. We believe the differences between 
the Balmer and Paschen dynamic spectra reflect their dif- 
fering optical depths; the Paschen lines are (relatively) less 
optically thick, and therefore come closer to tracing individ- 
ual clouds in the magnetosphere. 

Since we know the circumstellar material is bound in 
co-rotation, we can unambiguously map radial velocity onto 
the projected stellar surface, as indicated by the upper hor- 
izontal axis in Figs. [15] and |16| With this information, we 
can now measure the furthest extent of the magnetospheric 
emission from the centre of the star. We find that the emis- 
sion extends out to a maximum distance of 3.8 Req at phase 
0.39, 4.4: Req at phase 0.59, and 4.2 R^q at phase 0.93, as 
measured from the Ha dynamic spectrum. The middle of 
both redward features is found to occur at a distance of 
~1.7i?e<j, while the blueward feature is found to occur at a 
distance of 1.9 Req. 

We also estimated the density in the circumstell ar envi- 
ronment following the procedure of Stefl et al. ( 2003 1. Using 
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the BK3 model as a good approximation of the photospheric 
profiles, we estimated the Balmer decrements D34 and D54 
during phases where there is no apparent absorption in the 
dynamic spectra of Ha {(j> — 0.45, 0.72), from the resid- 
uals of the observed minus theoretical equivalent widths, 
corrected to absolute flux using /*(Ha)/f*(H/3) = 0.36 and 
/*(H7)/f*(H/3) = 1.44, as measured from the BK3 model. 
The Balmer decrements were then converted to densities 
using the theoretical computations of [Williams &i Shipman| 
1 1988 1, which were derived for an isothermal, pure hydrogen 



Table 3. Summary of stellar, magnetic, wind, and magneto- 
spheric properties of HR 5907 derived in this work. 



accretion disk of 10 000 K (the derived densities are weakly 
dependent on the temperature), which is optically thin in 
the continuum. 

From these measurements, we find that the logarithmic 
density ranges from 11.3 to 13.5 particles per cm^. There 
is a large discrepancy between the logarithmic densities de- 
rived using -D34, which give a value of 13.5 at both phases, 
compared to values computed from D54, (11.3, 11.6). We 
speculate that the difference may be due to the poor nor- 
malisation of the blue wing of the H/3 line and that the 
EW measurements of this line do not span the same veloc- 
ity range as the other Balmer lines, since part of the blue 
wing of this line is not covered over a single order in the 
ESPaDOnS spectra. Another possibility is that the magne- 
tospheric densities are very near the optically thick limit, or 
more likely straddling this limit, above which the measured 
Balmer decrements become independent of the density. 



8 DISCUSSION & CONCLUSIONS 

This paper reports extensive spectroscopic and mag- 
netic monitoring of the rapidly-rotating, early B-type star 
HR 5907, based on UVES and FORS data obtained at the 
VLT, ESPaDOnS data from the CFHT, and MOST pho- 
tometry. These observations and their analysis were under- 
taken within the context of the Magnetism in Massive Stars 
(MiMeS) Project. 

We combined newly measured MOST photometry with 
archival Hipparcos measurements to refine the photometric 
period of this star to ~0. 508276 d. A period search on the 
spectropolarimetric data confirms this period is also present 
in the equivalent width variations of Ha as well as the 
longitudinal magnetic field measurements. If we interpret 
these observations within the context of the ORM with a 
rigidly-rotating magnetosphere, then this period represents 
the rotation period of the star. A comparison of the Ha 
EW curve and photometric lightcurve suggests that there 
is ~0.05 phase offset between the minima. This phase off- 
set is not a result of a potentially incorrect period, as this 
would only affect the relative phasing between the Hippar- 
cos photometry and the current epoch of observations and 
this offset is still present between the MOST photometry 
and Ha EW curve, indicating that there is a real phase shift 
between the observed minima. As with the differences be- 
tween the Balmer and Paschen line profile variability, the 
phase shift likely arises from differences in optical depths of 
the continuum versus the Ha line. 

Using the optical spectroscopic data, archival lUE UV 
observations and the estimated rotation period, we were able 
to constrain the stellar physical parameters by comparing 



Stellar Parameters 

Stellar Mass (Mq) 
Teff (K) 
log(L*/L0) 
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~1.2 
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11.3-13.5 



|1997[ ), which properly reflect the oblateness of the star due 
to rapid rotation. The Teff we find (~17kK) is cooler than 
may otherwise be suggested by HR 5907's B2.5 spectral type, 
but can be understood if this star is actually a He-strong 
star. We believe this to be the case as the observed he- 
lium line depths are deeper than the modelled helium lines 



computed using solar abundances (as illustrated in Fig. 171 



This would make HR 5907 similar to other magnetic early B- 
type stars. We speculate that the helium over-abundance in 
HR 5907's photosphere is likely the reason that other recent 
studies ( Fremat et al.|[2"005 l that use helium lines as a pri- 
mary temperature indicator (even when taking into account 
rapid rotation) derive a much higher effective temperature 
(~21250K). We also investigated whether the discrepancy 
between the observed and modelled He lines was due to the 
assumption of LTE. Calculations using NLTE atmospheres 



the spectra to a grid of bruce/kylie models (Townsend 



(also displayed in Fig. 171 are not able to resolve the dis- 
crepancy. 

As a final test, we determined the effective temperature 
that would be necessary to approximately fit the observed 
He line depths. The derived value is about 2 kK hotter than 
our adopted temperature (i.e. 19 kK). Such a temperature 
results in the C 11 4267 A line being about 25 percent deeper, 
and the UV SED being about 20 percent too large. Carbon 
abundances as measured from optical lines are poorly stud- 
ied in Bp stars, but are typically found to be underabun- 
dant in cooler Ap stars (e.g. Roby & Lambert 19901. If 
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Figure 14. Phased equivalent width variations of hydrogen lines. 



Teff = 19 kK is adopted this would suggest that the carbon 
abundance would need to be reduced by 0.4 dex to fit the 
observed line depth, but this discrepancy will be somewhat 
mitigated if we use NLTE models, as the Cll 4267 A line 
is predicted to be weaker in NLTE, which would then in- 
crease the abundance necessary to fit the observed profile. 
However, adopting a NLTE model does not significantly af- 
fect the SED and the higher Tett would provide a poorer fit 
to the UV SED. A future self-consistent analysis that takes 
into account rapid rotation, NLTE effects, and non-solar He 
abundances is necessary to remove any ambiguity. 

From the low-resolution FORS and high-resolution ES- 
PaDOnS circular polarisation measurements, we find that 
the mean, surface averaged, longitudinal magnetic field val- 
ues are consistently negative. From the variation of these 
measurements we infer a dipolar magnetic field with a 
~16kG field strength at the pole and a magnetic axis 
nearly aligned with the rotation axis (/3 ~ 7°). However, 
direct modelling of the mean LSD profiles of the velocity- 
resolved Stokes V signatures suggests a lower dipole polar 
field strength of ~10kG. In any event, modelling of the cir- 
cular polarisation data shows large deviations from the sim- 
ple dipole model at several rotational phases. These impor- 
tant deviations may indicate a more complex topology to 
the magnetic field structure, similar to what is seen in some 
other magnetic B-type stars. However, the observed devia- 
tions could also be due to the strongly variable helium lines 
that dominate our mean LSD profiles. In fact, the lower Bt 
measurements around phase 0.5 are completely consistent 
with what is expected if helium is less abundant at the pole 
and therefore more weighted towards the equatorial regions 
where the helium lines are stronger. A self-consistent He- 
abundance and magnetic field map of the stellar surface, 
similar to HD 37776 (Kochukhov et al. 20111, is required 



to unambiguously determine the magnetic field topology of 
HR5907. 

Comparing the observed lightcurve with the photomet- 
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ric predictions of [Townsend (20081, using our inferred incli 



Figure 17. Comparison of observed helium lines (solid black; 
from 24 February 2010) with our best-fit solar abundance model 
(Tefi = 17kK, Mi, = 5.5M0, Req = 3.1 Rq) computed in LTE 
(dashed red) and NLTE (dotted green). As the Tog is indepen- 
dently constrained by the SED, this implies that HR5907 is a 
He-strong star. 



nation and magnetic obliquity, shows the two to be in good 
qualitative agreement. A lower inclination would not result 
in a perceivable dip in the lightcurve, unless the magnetic 
obliquity was about 30° or higher. On the other hand, if the 
inclination were higher, we would expect to see two dips 
in the lightcuve that would require a ~ld period (which 
would be inconsistent with the variations observed with the 
Be measurements) as opposed to the single dip that we do 
observe with our ~0.5d period. If the magnetic obliquity 
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Figure 15. Phased observational variations of the circumsteUar magnetosphere. Shown are the variations in the indicated Balmer lines 
after subtracting the best-fit BK3 photosphcric model (dashed red). The upper panels have a colour scheme chosen to show the full 
dynamic range of the emission variation, while the dynamic range in the lower panels are chosen to highlight the emission (green to red) 
features versus absorption (blue to black) features. Also shown are dashed vertical lines to indicate the rotational velocity (and radius) 
of HR5907 and dashed curves to highlight the motion of the emission features as discussed in the text. 
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Figure 16. Phased observational variations of the circumstellar magnetosphere. Shown are the variations in the indicated Paschen lines 
after subtracting the best-fit TLUSTY photospheric model (dashed red). The upper panels have a colour scheme chosen to show the full 
dynamic range of the emission variation, while the dynamic range in the lower panels are chosen to highlight the emission (green to red) 
features versus absorption (blue to black) features. Also shown are dashed vertical lines to indicate the rotational velocity (and radius) 
of HR5907 and dashed curves to highlight the motion of the emission features as discussed in the text. 
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was much larger, we would expect to see two dips in the 
lightcurve, unless the inclination was lower. 

As shown in Sect. [5] there appear to be significant chem- 
ical abundance variations across the stellar photosphere. He- 
lium appears to have the strongest contrast, with equivalent 
width variations on the order of 20 percent, which is actu- 
ally much smaller than observed in other rapidly-rotating 
He-strong stars such as a Ori E or HR7355. Carbon, sili- 
con, and nitrogen are also found to vary weakly. In other He- 
variable stars the locations of these abundance spots are cor- 
related with the magnetic field (e.g. [Bohlender &: Landstreet| 



1988) and usually axisymetric about the magnetic equator. 



Therefore the weak helium variability may well reflect the 
observed magnetic geometry - we are nearly always viewing 
the same orientation of the magnetic hemispheres and would 
expect only small rotational modulation since the magnetic 
axis is nearly aligned with the rotation axis. While the ab- 
solute helium variability might be weak compared to other 
rapidly-rotating He-strong stars, the helium surface distri- 
bution may vary enough to cause photometric fluctuations. 
The small dip that is evident in the photometric lightcurve 
(at about phase 0.65) also coincides with a small increase in 
helium absorption at this phase, as shown in Fig. [7] 

Using the stellar parameters as derived in this pa- 
per and wind parameters estimated from these stellar pa- 
rameters (included in Table [s] and obtained using formu 



lae from Castor, Abbott & Klein ( 1975 1 and Vink, de 
Koter & Lamers (2001 1) , we find that the wind mag- 



netic confinement parameter ( ud-Doula & Owocki 2002 1 
10' 



which indicates that the circum- 



stellar material is strongly confined out to large distances 
from the star (7?aif = i?*^* ~ 30 — 40 J?eq). The Kepler, or 
co-rotation, radius ( |ud-Doula, Owocki &: Townscnd 2008] ) 



3/2 0;-^/=* 7?^ 



■pole. 



2.0 Req, implying that beyond this 



radius, the material is centrifugally supported. We find that 
the locations of the emission peaks in Balmer and Paschen 
lines agree well with this value. We also computed the pre- 



dicted spin-down time using Eq. (25) of ud-Doula, Owocki 
|fc Townsend| p009t , Tspin ~ 8Myr, which is considerably 
longer than the potential age of HR 5907. 

In summary, this study reports the detection of a large- 
scale, organised magnetic field with a polar surface intensity 
of 10-16 kG in HR 5907. Interpreted within the context of the 
ORM, we find photometric, Hq and longitudinal magnetic 
field variations consistent with a 0.508276 d rotational pe- 
riod, making this the shortest period, non-degenerate, mag- 
netic massive star known to date. We find evidence of surface 
abundance variations in helium, carbon, silicon, and nitro- 
gen. 

This is one of only a few known rapidly-rotating mag- 
netic massive stars that show strong emission variations due 
to a magnetosphere. The only other massive star with a 
comparable rotation period is HR7355, which is believed to 
have a magnetic geometry more similar to a Ori E (a mag- 
netic obliquity closer to 90°), which is very different from 
this star. HR5907 is therefore an ideal target for compari- 
son with the predictions of rigidly rotating magnetosphere 
models ( |Townsend fc Owocki||2005| |Townsend||2008| ) and a 
great testbed for studying the effects of the magnetic field 
orientation on angular momentum loss and magnetic spin- 
down. 
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